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ABSTRACT 

We show that primordial nucleosynthesis in baryon inhomogeneous big-bang mod- 
els can lead to significant heavy-element production while still satisfying all of the light- 
element abundance constraints including the low lithium abundance observed in popula- 
tion II stars. The parameters which admit this solution arise naturally from the process of 
neutrino induced inflation of baryon inhomogeneities prior to the epoch of nucleosynthe- 
sis. These solutions entail a small fraction of baryons (^2%) in very high density regions 
with local baryon to photon ratio r\ h ~ 10 -4 , while most baryons are at a baryon-to- 
photon ratio which optimizes the agreement with light-element abundances. This model 
would imply a unique signature of baryon inhomogeneities in the early universe, evidenced 
by the existence of primordial material containing heavy-element products of proton and 
alpha-burning reactions with an abundance of [Z] ~ —6 to —4. 

Subject headings: early universe - abundances, nuclear reactions, nucleosynthesis 



1 



1. Introduction 

Recently, there has been considerable interest in the possibility that primordial nu- 
cleosynthesis may have occurred in an environment in which the baryons were distributed 
inhomogeneously (Alcock, Fuller, & Mathews 1987; Applegate, Hogan, & Scherrer 1987; 
1988; Fuller, Mathews, & Alcock 1988; Kurki-Suonio et al. 1988; 1990; Malaney & Fowler 
1988; Boyd & Kajino 1989; Terasawa & Sato 1989abc; 1990; Kajino & Boyd 1990; Kurki- 
Suonio & Matzner 1989; 1990; Mathews et al. 1990, 1993; Kawano et al. 1992; Jedamzik, 
Fuller & Mathews 1993; Thomas et al. 1993). Such studies were at first motivated by 
speculation (Witten 1984; Applegate & Hogan 1985) that a first-order quark-hadron phase 
transition at (T ~ 100 MeV) could produce baryon inhomogeneities as baryon number was 
trapped within bubbles of shrinking quark-gluon plasma. Even if this transition is not first 
order as some lattice-gauge theory calculations now seem to imply (e.g. Brown et al. 1990) 
such fluctuations might also have been produced by a number of other processes operating 
in the early universe such as late baryogenesis during the electroweak phase transition 
(Fuller et al. 1993), a TeV-scale Z(3) QCD phase transition (Ignatius, Kajantie, & Rum- 
mukainen 1992), kaon condensation after the QCD transition (Nelson 1990), or magnetic 
fields from superconducting cosmic strings (Malaney & Butler 1989). Other mechanisms 
have also been proposed (see Malaney & Mathews (1993) for a recent review). 

The first calculations (Alcock et al. 1987; Applegate et al. 1987; 1988; Malaney & 
Fowler 1988; Fuller et al. 1988) of the effects of such inhomogeneities were particularly 
interesting since they seemed to admit a much larger average universal baryon-to-photon 
ratio rj than that allowed in the standard homogeneous big bang model (Wagoner, Fowler, 
& Hoyle 1967; Schramm & Wagoner 1977, Yang et al. 1984; Boesgaard & Steigman 
1985; Olive et al. 1990; Walker et al. 1991; Smith, Kawano & Malaney 1993) while 
still satisfying the light-element abundance constraints. Subsequent calculations (Kurki- 
Suonio et al. 1990; Malaney & Fowler 1988; Terasawa & Sato 1989abc; 1990; Kajino & 
Boyd 1990; Kurki-Suonio & Matzner 1989; 1990; Mathews et al. 1990, 1993; Kawano et 
al. 1992; Jedamzik et al. 1993; Thomas et al. 1993) have shown, however, that when 
baryon diffusion is properly coupled with the dynamics of the system, the allowed values 
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of rj which are consistent with the light-element abundances are not significantly different 
than those of the standard big bang. 

Even with this constraint on i], however, an interest in baryon inhomogeneities has 
remained due to the possibility that such fluctuations might have produced an observable 
signature in the abundances of heavier elements such as beryllium and boron (Kajino 
k Boyd 1990; Malaney k Fowler 1989; Terasawa k Sato 1990; Kawano et al. 1991), 
intermediate mass elements (Kajino, Mathews k Fuller 1990) or heavy elements (Applegate 
et al. 1988; Applegate et al. 1993). Such possible signatures are also constrained, however, 
by the light-element abundances (e.g. Terasawa k Sato 1990). In particular, the population 
II lithium abundance, Li/H ~ 10~ 10 , constrains the possible abundances of synthesized 
heavier nuclei to be quite small (e.g. Alcock et al. 1990; Terasawa k Sato 1990). 

The purpose of this paper, however, is to point out a previously unexplored region 
of the parameter space for baryon-inhomogeneous big bang models in which substantial 
production of heavy-elements is possible without violating any of the light-element abun- 
dance constraints. This solution arises naturally from one of the most likely scenarios for 
the production of baryon inhomogeneities in the first place. In fact, we discovered this 
solution not by looking for ways to optimize heavy-element production, but by exploring 
the natural evolution of large baryon number-density fluctuations. If such fluctuations 
are created in an early epoch then they will be damped subsequently by neutrino induced 
heating and expansion to a characteristic baryon-to-photon ratio ~ 10~ 4 . (Jedamzik k 
Fuller 1993; Jedamzik et al. 1993). 

2. Baryon Inhomogeneities and Neutrino Inflation 

Let us begin with a discussion of how baryon inhomogeneities might be generated 
in the early universe and how they are most likely to appear. The most likely scenario 
probably involves a first-order phase transition. It is not clear whether that transition is 
associated with the QCD epoch or an earlier time. Nevertheless, any process which couples 
to baryons in a spatially-dependent way can produce baryon inhomogeneities (Malaney k 
Mathews 1993). For example, in the QCD transition, inhomogeneities would most likely be 
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produced through a combination of the limited permeability of the phase boundary and the 
higher thermodynamic solubility of baryon number in the high-temperature phase (Witten 
1984; Fuller et al 1988; Kurki-Suonio 1988). The same might be true if inhomogeneities 
are produced by a kaon condensate (Nelson 1990) after the QCD transition. 

If baryogenesis is associated with a first-order phase transition at the electroweak 
epoch (Kuzmin, Rubakov, & Shaposhnikov 1985; Cohen, Kaplan & Nelson 1990; Dine et 
al. 1991; Turok & Zadrozny 1991; McLerran et al. 1992) then it may well be that the net 
baryon number is generated in a inhomogeneous fashion (Fuller et al. 1993). The pres- 
ence of cosmic strings might also induce baryon inhomogeneities through electromagnetic 
(Malaney & Butler 1989) or gravitational interactions. 

It is a common misconception that the most natural amplitude for baryon-number 
fluctuations is just that given by the thermodynamic ratio of equilibrium baryon densities 
in the two phases of a QCD transition. This as a number ~ 100 for the QCD transition 
(Alcock et al. 1987). This value for the fluctuation amplitude, however, is unlikely. It 
would occur only if complete equilibrium were maintained in both phases till near the end 
of the phase transition followed by a sudden complete drop from equilibrium. This would 
require efficient mixing of baryon number in both phases and efficient transport of baryon 
number across the phase boundary till just near the end of the phase transition. Near the 
end of the phase transition, the transport of baryon number across the phase boundary 
can not be efficient enough to establish chemical equilibrium. This is because the velocity 
of the phase boundary continuously increases. Efficient mixing of baryon number within 
the two phases will also be unlikely due to the short mean-free path for baryons for T > 
50 MeV. Furthermore if baryon mixing were this efficient throughout the phase transition, 
it would likely be sufficiently efficient after the phase transition to homogenize any baryon 
inhomogeneities long before the epoch of nucleosynthesis. 

Mechanisms which take into account the efficiency for baryon transport, in partic- 
ular, can lead to significant baryon inhomogeneities, and even in the most conservative 
conditions (e.g. Kurki-Suonio 1988), can produce a central baryon density far in excess 
of the thermodynamic equilibrium density (Fuller et al. 1988). If the baryon fluctuations 
are produced by such a mechanism, there is almost no limit as to how high the central 
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baryon-to-photon ratio can become. The possibility of such a high baryon density in fact 
motivated Witten (1984) to suggest this mechanism as a means to form stable quark- 
matter agglomerates during the QCD phase transition if such matter were the ground 
state of baryonic matter at high density. 

In any event, it is certainly possible that baryon inhomogeneities produced by any of 
the above mechanisms could reach a central baryon-to-photon ratio of 77 > 10~ 4 . This 
being the case, it is important to note the effect of neutrinos on any such large-amplitude 
baryon fluctuation (Hogan 1978; Hogan 1988; Heckler & Hogan 1993; Jedamzik & Fuller 
1993). These over dense regions will have a higher pressure contribution from baryons than 
their surroundings. Hence, in order to maintain pressure equilibrium with the surrounding 
photon bath, any region of baryon overdensity will be at a slightly lower temperature than 
the baryon depleted regions. However, since neutrinos have a longer mean free path than 
photons they will penetrate into and heat the high-baryon-density regions causing them 
to expand and diminish the excess baryon density. 

Several recent papers have explored this process in detail (Heckler & Hogan 1993; 
Jedamzik & Fuller 1993; Jedamzik et al. 1993). Jedamzik & Fuller (1993) found that 
this neutrino-induced damping of inhomogeneities is more pronounced as the density in- 
creases, such that any fluctuation with baryon-to-photon ratio in the high-density regions, 
rj h > 10 -4 , is inflated to rj ~ 10 -4 . This is independent of its initial baryon density and 
holds true over a wide range of initial fluctuation radii. Fluctuations with rj less than this 
are essentially unchanged. 

This is illustrated in figure 1 which shows the final baryon-to-photon ratio ?]/ after 
neutrino inflation as a function of initial comoving radius of a fluctuation. Lines are drawn 
corresponding to different initial baryon-to-photon ratio, r/j. The degree to which the 
fluctuations are damped is related to the time scale during which the neutrino mean free 
path is of order the radius of the fluctuation. For larger fluctuation radii, this time scale 
becomes shorter due to the rapid increase in neutrino mean free path as the temperature 
decreases. Hence, if the initial radius of the fluctuation is too large, neutrinos can not 
sufficiently damp the fluctuations to bring them to the limiting value of t] ~ 10~ 4 . 
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The effects of neutrino inflation on very-large amplitude baryon fluctuations with ini- 
tial baryon-to-photon ratios of order unity (rji ~ 1) are not shown in figure f . Fluctuations 
at such high baryon densities could correspond to stable or metastable low-entropy rem- 
nants like strange-quark matter nuggets (Witten 1984). In this case heating by neutrinos 
at temperatures of T 50 MeV will cause the evaporation of such nuggets rather than 
inflation. Only the very largest nuggets will survive complete evaporation (Alcock & Farhi 
1985; Madsen, Heiselberg & Riisager 1986). Evaporated nuggets can leave behind a fluctu- 
ation with a baryon-to-photon ratio significantly larger than rj ~ 10 -4 even after neutrino 
inflation. The exact baryon distribution after evaporation will depend on the initial baryon 
number of the nuggets as well as on bulk properties of quark matter and strong interaction 
physics at the surface of the nuggets. 

In any case, fluctuations with small initial radius and initial baryon-to-photon ratio 
less than unity will be inflated to a natural amplitude of rj ~ 10 -4 by the time of primor- 
dial nucleosynthesis. The average baryon-to-photon value in the baryon-depleted regions 
rj 1 (or equivalently the total average baryon-to-photon ratio, rj) and the total fraction of 
baryons in the high density regions, can then be fixed by the model for forming the 
fluctuations or treated as free parameters to be fixed by the constraints from light-element 
nucleosynthesis. As noted above, models for the formation of inhomogeneities are consis- 
tent with small values of We will now show that a small value for is also desired for 
agreement with light-element nucleosynthesis. 

3. Results 

Fixing r\ at 10 -4 in the high-density regions can avoid the overproduction of 7 Li. This 
was apparent already in the standard big bang calculations of Wagoner et al. (1967). 
This value for 1] is, however, much larger than has been considered in most standard or 
inhomogeneous big bang models since that early work. We can illustrate the effects of a 
high rj h value by first considering the simplest possible baryon inhomogeneous model, that 
in which the fluctuations are so widely separated that no significant baryon diffusion occurs 
before or during the epoch of primordial nucleosynthesis. In this limit the high-density 
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and low-density regions can be treated (Wagoner et al. 1967; Wagoner 1973; Schramm 
& Wagoner 1977; Yang et al. 1984) as separate standard homogeneous big bang models 
with different local values for r\. The nucleosynthesis yields from these regions can then 
be averaged after the epoch of nucleosynthesis. (We will also consider below the effects of 
baryon diffusion on small-scale fluctuations.) 

The results of standard big bang nucleosynthesis (with three light neutrino species) 
at such high densities are summarized in figure 2 where we have extended the big-bang 
nuclear reaction network through mass A = 28. The point which has been overlooked in 
recent studies of baryon inhomogeneous models is that, for r\ greater than a few times 
10~ 5 , 7 Li is destroyed by the 7 Li(a,7) 11 B reaction. The 11 B thus formed is then rapidly 
consumed by a sequence of proton and alpha captures up to heavy nuclei reminiscent of 
an rp-process (Wallace & Woosley 1981). In fact, for 1] > 10~ 4 , almost all of the material 
with atomic mass heavier than 4 He is processed to nuclei beyond the end of our network 
at A = 28. Thus, we avoid the overproduction of 7 Li usually associated with averaging 
over i] values near the lithium minimum at r\ ~ 2 x 10~ 10 (e.g. Alcock et al 1987). 

It is, of course, still necessary to produce some deuterium, and in so doing, some 
3 He and 7 Li will also be produced. There is also the problem that the 4 He abundance 
is significantly overproduced, i.e. Y p ~ 0.36 for rj ~ 10~ 4 . The simultaneous solution to 
these constraints fixes allowed values of rf and the fraction of baryons in the high-density 
region, /& = f v (v h /v) where f v is the fraction of total volume occupied by the high 
density regions, and r] is the average baryon-to-photon ratio. 

We adopt constraints on observed light-element abundances from Walker et al. (1991) 
and Smith et al. (1993) such that in this simple model with no diffusion, the set of 
constraint equations become, 



0.22 < hYh + il-fjYt < 0.24 



(la) 



D/H 




-5 



(16) 



[D+ 3 He]/H ~ (1 - f b )([D+ 3 He]/H) 1 < 1.0 x 10" 4 



(lc) 



0.8 x 10 



-10 



< Li/H 



(1 - f b )(Li/Hy < 2.3 x 10 



-10 



(Id) 
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where Xi is the mass fraction of species, i, and the superscripts h and I refer to quantities 
computed in the high and low baryon density regions, respectively. The approximate 
relations in (lb)-(ld) are valid for small fb and insignificant D, 3 He or 7 Li production in 
the high baryon density regions. 

Figure 3 shows allowed regions of the parameter space of rf , (in units of 10~ 10 ) as 
a function of fb- Clearly there is a region of fb and rf for which all of the light-element 
abundance constraints can be satisfied. This region is mostly defined by the upper limits 
to helium and [D + 3 He]/H. This fi sure does not include the extension of the allowed 
parameter space due to nuclear-reaction uncertainties, as is usually done in standard big- 
bang studies (e.g. Walker et al. 1991; Smith et al. 1993). For example, if we had allowed 
for uncertainties in helium and deuterium production, the range of allowed rj[ would 
expand to 2.8 < rj[ < 4.0 for small fb, This is just the usually- quoted standard big- 
bang limit for 3 light neutrino species. The most stringent constraint on fb is from the 
the upper limit to the primordial helium abundance which necessitates that at most only 
~ 2% of the baryons can be in high density regions. If the uncertainty in the neutron half 
life is also included, this upper limit increases to as high as /t ^ 0.1. 

It is easy to see from an inspection of Figure 2 why this region of the parameter space 
is allowed. By placing almost all of the baryons (^98%) in a low density region for which 
rj is close to the optimum standard big bang value, a good fit to all of the light-element 
abundances is guaranteed. 

So we see that in this large-separation limit for inhomogeneous models, all of the 
light-element constraints can be satisfied. Even though this constraint requires that only 
a small fraction of the baryons be in high-density regions, we still find a potentially in- 
teresting (though not necessarily practically observable) signature from inhomogeneous 
nucleosynthesis. For r\ h ~ 10~ 4 and fb ~ 0.01, the mass fraction Z of heavy (A > 28) or 
CNO nuclei is Z ~ 3 x 10~ 9 ([Z] ~ —6.5) after averaging with the low-density regions. 

The mass fraction in nuclei beyond the end of our network at A = 28 could reach 
as high as [Z] ~ —4 for rf 1 ~ 10 ~ 3 and ~ 0.01. A baryon-to-photon ratio as high 
as rj ~ 10~ 3 after the epoch of neutrino inflation requires an initially large-amplitude 
fluctuation of large radius, or equivalently a large baryonic mass M b ~5x 10 _11 M Q (cf. 
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Fig.l). Another possibility to obtain a local baryon-to-photon ratio during primordial nu- 
cleosynthesis in excess of rf 1 ~ 10~ 4 is from the evaporation of low-entropy nuggets. The 
effects of strange-quark matter nuggets on primordial nucleosynthesis have been considered 
(Schaeffer, Delbourgo- Salvador & Audouze 1985; Madsen & Riisager 1985), but the evapo- 
ration of nuggets before primordial nucleosynthesis and the possibility of significant heavy 
element production during primordial nucleosynthesis has not been addressed before. 

It is interesting to note that a small fraction of baryons initially within large-mass 
fluctuations or low-entropy nuggets might produce a mass fraction in heavy elements during 
primordial nucleosynthesis, in excess of [Z] ~ —4. Thus, observation (e.g. Beers, Preston, 
& Shectman 1992) of very low-metallicity objects {[Z]< — 4) can potentially constrain 
the existence of large-mass fluctuations or low-entropy remnants in the early universe. 
The existence of a floor in primordial heavy-element abundances could even provide a 
nucleosynthesis signature of the physical properties of the fluctuations. 

In any case, fluctuations of moderate amplitude and size will be damped by neu- 
trino inflation to a baryon-to-photon ratio of r\ h ~ 10~ 4 . Hence, [Z] ~ —6 is probably 
a good characteristic abundance for the sum of heavy-element abundances produced by 
inhomogeneous nucleosynthesis. 

In this possible abundance diagnostic signature includes C/O and N/O ratios 

(Kajino et al. 1990). A high C/O or N/O ratio is contrary to what one might expect from 
a first generation of massive stars, but does appear for big bang models with 10~ 5 < i] h < 
10~ 4 . For rj h > 10~ 4 , one should see an anomalously low [O/Fe] ratio, assuming that the 
A > 28 nuclei are burned to Fe. We realize, of course, that detecting abundances this low 
is extremely difficult and probably not possible with existing techniques. Nevertheless, it 
may eventually be possible to accurately measure abundances this low, in which case the 
presence or absence of baryon inhomogeneities might be explorable. 

3.1 Effects of Baryon Diffusion 

The results shown on figure 3, of course, correspond to only one possible model in 
which the separation between fluctuations is so large that baryon diffusion is insignificant. 
We have also explored models in which neutron diffusion has been included but for a 
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reaction network limited to A < 12. We have considered models in which the proper (at 
T = 100 MeV) separation distance was varied between the QCD horizon (10 4 m) and a 
small fraction thereof (0.1 m) for fluctuations with fixed r\ h and rj . Alternatively, we have 
considered effects of baryon diffusion as a function of baryon mass in the fluctuation, Mb, 
and fraction ft, of the baryons in the high density regions for fixed rf . 

Conventionally a regular lattice of fluctuations is described by four variables, the mean 
separation between adjacent fluctuations, r, as well as three more variables. These can 
be chosen as the volume fraction in the high-density region, /„, the ratio of densities in 
high-density region and low-density region, R , and the average baryon-to-photon ratio, r] 
(Fuller et al. 1988). In the present analysis of high-amplitude fluctuations, which include 
only a small fraction of the total baryons ( ft, « 1), it is convenient to characterize the 
fluctuation by baryon mass instead of the mean separation distance. Varying fb at fixed 
Mb is equivalent to varying the separation distance. Varying Mb at fixed fb corresponds 
to a self-similar variation of fluctuation radius and separation distance. 

A conversion between the mean separation and baryon mass is easily obtained, from 



where r is the proper mean separation between fluctuations (at T = 100 MeV), and M b 
is the baryon mass within a high-density fluctuation. In equation (2) and what follows 
it is straightforward to convert from proper length at To = 100 MeV to proper length at 
temperature T. To a good approximation, one can do this by simply assuming that the 
scale factor a oc Tq/T for all temperatures T > 0.1 MeV. The baryon mass is related to 
the initial (final) fluctuation baryon-to-photon ratio, r\i (?7/), and fluctuation radius, Ri 
(Rf), as follows 



where the proper (T = 100 MeV) fluctuation radius R is in meters. Note that the prod- 
uct of baryon-to-photon ratio times the cube of the fluctuation radius is not changed by 
neutrino inflation because baryon number within the fluctuation is conserved. 





M^4x lO" 13 M ?7f Rf = 4 x 10~ 13 MqT] 1 } R 3 f 



(3) 
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With the help of equation 3 and figure 1 the reader is able to determine the mass of 
a fluctuation, its baryon-to-photon ratio after neutrino inflation, rjf , and radius, Rf, 
assuming rji and Ri. 

In order to identify the mass scales for which baryon diffusion affects fluctuations 
before and during primordial nucleosynthesis, the baryon diffusion length, d, should be 
compared to the radius of the fluctuation after neutrino inflation, Rf. Baryon diffusion 
in high-density regions is limited by neutron-proton scattering (Applegate et al. 1987, 
Jedamzik & Fuller 1993). The proper (referred to Tq = 100 MeV) baryon diffusion length 
at temperature T, i.e. the average distance a neutron diffuses between high temperatures 
and temperature T , is given by 



Clearly, the baryon diffusion length is smaller for a larger value of the baryon-to-photon 
ratio in the fluctuation. For a fluctuation of baryon mass, M b <1O~ 21 M , the baryon 
diffusion length at T ~ 500 keV is comparable to the fluctuation size, and the fluctuations, 
therefore, vanish by diffusion before nucleosynthesis begins. 

For mass scales Mf, > 10 -21 M©, for which the fluctuations can survive to the epoch 
of nucleosynthesis, the abundance yields from models which include baryon diffusion are 
not that much different from figure 2. This is because in high-density regions with r\ h ~ 
10~ 4 , nuclear reactions can begin to synthesize light elements at temperatures as high 
as T 300 — 400 keV. Thus, in the high-density regions free neutrons are incorporated 
into helium soon after the freezeout of weak reactions. As the weak reactions freeze out, 
neutrons and protons are no longer rapidly converted from one into another. Hence, they 
obtain separate identities. 

Because the neutrons are absorbed so quickly after weak freezeout, there is not much 
time for neutron diffusion to produce large variations in the neutron-to-proton ratio. Even 
though some free neutrons in the low-density regions can diffuse back into the high-density 
regions (Malaney & Fowler 1988), this back diffusion does not significantly affect the light 
element abundances when << 1. This is because light elements are produced in low- 
density regions. This influx of neutrons into the high-density regions may somewhat affect 
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the yields of heavy elements which are produced there. For the most part, however, the 
neutron diffusion length in the high-baryon-density regions is short and the final abundance 
yields should not differ much from that obtained by simply averaging over regions with 
different r\ and fixed n/p ratio as in figure 2. 

The results of numerical studies including the effects of baryon diffusion are presented 
in figures 4 and 5. These figures illustrate that fluctuations more widely separated than a 
proper distance of r = 100 m (at T = 100 MeV) or of a mass scale > 10~ 18 M Q produce 
nucleosynthesis yields which are indistinguishable from those of figure 3 (Mathews et al. 
1990). Similarly, models in which the proper separation distance (at T = 100 MeV) 
is r<5m (M b < 1O _21 M ) are indistinguishable from a homogeneous standard big bang 
model with the same average r\. The reason is that baryon diffusion prior to nucleosynthesis 
damps out the inhomogeneities (Mathews et al. 1990). The maximum effects on the 
nucleosynthesis yields occur for proper separation distances (at T = 100 MeV) between 
10 and 60 m where the separation roughly equals the neutron diffusion length during 
primordial nucleosynthesis. (Applegate et al. 1988; Mathews et al. 1990; Terasawa & Sato 
1990). 

Figure 4 identifies the region in the Mb (total baryon mass of the fluctuation) versus 
fb plane in which heavy-element nucleosynthesis may imply a signature of high-amplitude 
fluctuations. To produce this figure we used the inhomogeneous nucleosynthesis code of 
Jedamzik et al. (1993) which includes hydrodynamic and radiation transport effects. We 
evolved baryon-number fluctuations between T = 100 MeV and the end of primordial 
nucleosynthesis at T ~ 5 keV. The effects of neutrino inflation at high temperatures are 
included in the numerical treatment. 

The total baryon mass within the horizon for the electroweak transition (M EW ) and 
the QCD transition (Mqcd) is also identified on figure 4. For Mb >10~ 18 (of order 
the electroweak baryon mass scale) and fb ^0.02 heavy elements are formed without 
overproducing lithium or helium. For less massive fluctuations, M b ~ 10~ 21 , it is also 
possible to satisfy the lithium and helium constraints, but in this case baryon diffusion 
erases the inhomogeneities before the onset of primordial nucleosynthesis so that the results 
are indistinguishable from the homogeneous case. 
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The effects of diffusion on the yields of light abundances is not very large even at 
optimum separation distances or masses. This is because the high-density regions are 
constrained to contain only a small fraction of the total baryons. The main contribution 
to the light-element abundances stems from the low-density regions which are not much 
affected by diffusion from the high-density regions. 

The effects of baryon diffusion on the various light-element abundances are briefly 
summarized as follows: 

Y P 

Independent of rj h , the helium mass fraction decreases by AY p ~ 0.0005 as the 
separation is decreased from 100 to 6 m. [For the remainder of this section all quoted 
length scales are proper lengths measured at T = 100 MeV.] As the separation is 
further decreased from 6 to 3 m, the helium mass fraction decreases by AY P ~ 0.002. 
This change is small compared to the uncertainty in the upper limit to Y p . 

D/H, (D+ 3 He)/H 

Because these light elements are made almost exclusively in the low-density regions 
which are largely uneffected by baryon diffusion their only dependence upon the dif- 
fusion length is to decrease by about 4% as the diffusion length is varied between 
50 and 5 m. Again, this change is insignificant compared to the uncertainties in the 
abundance constraints themselves. 

7 Li 

Lithium is the light-element most significantly affected by neutron diffusion. This is 
largely because of 7 Be production at the fluctuation boundary (Mathews et al. 1990). 
This causes the lithium abundance to increase to a maximum of 7 Li/H ~ 2 x 10~ 10 
for a separation of 10 m. For distances less than 6 m or greater than 100 m, the 
7 Li abundance remains at the value on figure 1 corresponding to whatever rf value is 
chosen. 

9 Be, 11 B 
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Since there has been some discussion of the production of these light nuclei in recent 
literature, a mention of their production in the present models should be made here. 
There can be no production of Be in these models, and almost no 11 B production 
[X( 11 B) <10 -18 ] even with an optimized separation distance (which, like lithium, can 
only increase boron by less than a factor of 2). Therefore, the abundances of these 
elements can not be a diagnostic for the kind of inhomogeneities considered here. 

A>12 

By far the yield which is most sensitive to the diffusion efficiency is that of the heavy- 
element abundances. This sensitivity is illustrated on figure 5. For separation dis- 
tances less than ~ 6 m the yields are no different than a standard homogeneous big 
bang with the same average value of rj. In this case most of the yield is in 12 C at an in- 
significant abundance relative to hydrogen of less than 10~ 14 (see figure 2). However, 
as the diffusion length is increased to ~ 100 m, the heavy-element abundances quickly 
increases by 5 orders of magnitude due largely to alpha and proton captures in the 
high-density regions. For r^^lO -4 , heavy elements are mostly produced in the form 
of iron and an anomalously low [O/Fe] would be a diagnostic abundance signature. 
For lO^^ry^lO -4 , anomalously high [C/O] and [N/O] ratios are expected (Kajino 
et al. 1990). Any observation suggesting the presence of primordial heavy elements 
with abundances as large as [Z] > - 6 could be a strong verification of the presence of 
baryon inhomogeneities in the early universe. 

4. Conclusions 

We have shown that there exists a plausible region of the parameter space for baryon- 
inhomogeneous primordial nucleosynthesis models in which all of the light-element abun- 
dance constraints are satisfied, yet which also leads to the unique production of potentially 
observable abundances of heavy nuclei. Although the abundances of these heavy nuclei are 
small, they can be as much as 6 to 8 orders of magnitude larger than in the homogeneous 
big bang. The region of the parameter space in which heavy-element abundances are gen- 
erated is naturally motivated by models for the formation of baryon inhomogeneities, for 
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example during a first-order phase transition, and by numerical simulations of the effect 
of neutrino-induced heating and expansion of baryon inhomogeneities prior to the epoch 
of primordial nucleosynthesis. Hence, this region of the parameter space should be taken 
seriously. 

The conditions most likely to result from the expansion of a high baryon density region 
by neutrino inflation which also satisfy the light-element abundance constraints, imply that 
a mass fraction as large as [Z] ~ —4 to —6 could be in the form of nuclei with A > 28 
and/or CNO nuclei after primordial nucleosynthesis. Thus, a search for a lower limit to 
the abundance ratios of CNO or heavy nuclei to the light primordial abundances could 
some day be a definitive indicator of the presence or absence of large-amplitude baryon 
inhomogeneities in the early universe. 
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Figure Captions 

Figure 1 Final baryon-to-photon ratio, rjf of neutrino-inflated fluctuations as a function of 
initial fluctuation radius (comoving at 100 MeV). Results are shown for four different 
initial values of the baryon-to-photon ratio, rji. 

Figure 2 Nucleosynthesis yields as a function of baryon-to-photon ratio, r\ for a reaction network 
extending from light nuclei to A = 28. These yields are for locally homogeneous regions 
sufficiently separated that baryon diffusion is insignificant. Note that for r\ > 10~ 4 
the yields of nuclei heavier than helium are almost entirely in massive nuclei with 
essentially no production of other light elements. 

Figure 3. Light-element abundance constraints on the baryon-to-photon ratio in the low density 
regions, r\ l ', as a function of the fraction, f b , of baryons in high-density regions. Note 
that both quantities, rf and f b , are plotted on a logarithmic scale. For this model the 
high baryon density region is fixed by neutrino induced inflation of the fluctuations 
to be?/ 1 = 10" 4 . 

Figure 4. Heavy-element nucleosynthesis signatures of high-amplitude fluctuations [r\i > 10~ 4 , 
rf = 2.8 x 10~ 10 ). Identified are regions of the total baryon mass M b (in units of M ) 
contained in a single fluctuation versus the fraction of baryons in these high-density 
fluctuations, f b . Also identified on this figure are the total baryon masses within the 
horizon for the electroweak transition (Mew) an d the QCD transition (Mqcd)- For 
M b > 10~ 18 (roughly the electroweak baryon mass scale) and f b ^0.02, heavy elements 
are formed without overproducing lithium or helium. For less massive fluctuations, 
M b < 10~ 21 , it is also possible to satisfy the lithium and helium constraints, but in this 
case baryon diffusion erases the inhomogeneities before the onset of primordial nucle- 
osynthesis so that the results are indistinguishable from the homogeneous standard 
big bang and no heavy-element signature emerges. 
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Figure 5. An example of the effects of baryon diffusion on heavy-element abundances (A > 
12). Shown are averaged yields for heavy elements as a function of proper separation 
distance in m at T = 100 MeV. For this figure we have arbitrarily set /& = 0.0001, 
rj = 2.8 x 10 -10 , and rj h = 1.3 x 10~ 5 . At separations less than 6 m, the results are 
indistinguishable from a standard homogeneous big bang with r\ ~ 2.8 x 10~ 10 because 
baryon diffusion before nucleosynthesis homogenizes the fluctuation. For separation 
distances greater than ~ 100 m, the results are essentially identical to figure 2. A 
transition between the two limits occurs when the neutron diffusion length during 
nucleosynthesis is comparable to the separation distance. 
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